On how a joint interaction of two innocent partners (smooth advection &c linear 

damping) produces a strong intermittency 



M. Chertkov 

Physics Department, Princeton University, Princeton, NJ 08544 
(February 5, 2008) 

Forced advection of passive scalar by a smooth d-dimensional incompressible velocity in the 
presence of a linear damping is studied. Acting separately advection & dumping do not lead to an 
essential intermittency of the steady scalar statistics, while being mixed together produce a very 
strong non-Gaussianity in the convective range: q-th (positive) moment of the absolute value of 
scalar difference, (\0(t;r) -6(t;0)\ q ) is proportional to r 5 ", £ 9 = y/d 2 /A + adq/ [(d - 1)D] - d/2, 
where a/D measures the rate of the damping in the units of the stretching rate. Probability density 
function (PDF) of the scalar difference is also found. 
Advection of passive scalar 8(t; r) by incompressible velocity field is a classical problem in turbulence theory. The 
problem attracts a lot of recent attention for remarkable combination of a rich physics and nontrivial rigorous results 
derived. R.H. Kraichnan pioneered the rigorous study of the problem inventing the temporal short-correlated but 
spatially non-smooth model of velocity for which the simultaneous pair correlation function of the scalar was found 
ty. However, the question of possible anomalous behavior of higher order (n > 1) structure functions S^^r) = 

\(9(t; r) — 9(t; 0)) 2 ™^} ~ r^ 2n was posed only 25 years later Next, the anomalous scaling A 2 „ = n£ 2 — £an> 
describing the law of the algebraic growth with L jr (where L is the scale of the scalar pumping) of the dimensionless 
ratio S2n(r)/ [^(r)]™, was shown to exist generically The anomalous exponent was calculated perturbatively 

in expansions about three non- anomalous ( A 2 „ = 0) limits, of large space dimensionality d ||||, of extremely 
non-smooth and almost smooth || velocities respectively. A strong anomalous scaling (saturation of £ 2 „ to a 
constant) was found for the Kraichnan model at the largest n by a steepest descent formalism ||. Although the 
£f) 1 restricted asymptotic information about anomalous exponent in the model is available a future possibility to establish 
a whole rigorous dependence of £ 2n on n, d and degree of velocity non-smoothness seems very unlikely (in a sense, the 
recent Lagrangian numerics H on the £4 as a function of (2 at d = 3 compensates the lack of rigorous information). 
In the present letter I discuss yet another passive scalar model with nontrivial anomalous behavior, £ 2 „ < n£; 2 , which 
£h ■ is possible to resolve explicitly for all the values of the governing parameters. The model describes generalization 
of the smooth (Batchelor) limit of the Kraichnan model on the case of a linear damping of the scalar. The pure 
Batchelor model (no damping), studied detaily in ^j, shows non-anomalous (logarithmic and almost Gaussian 
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for the Gaussian form of the pumping) behavior. The limit of a huge linear damping (neglect advection) is also non- 
anomalous. However, as it will be shown below (see (^,0)), the anomalous scaling does exist generically, and appears 
to be a nontrivial function of n, d and a parameter standing for the damping-to-convection ratio. The model describes 
forced advection of a scalar pollutant in the viscous-convective range (the viscous-to-diffusivity ratio is supposed to 

. £L ' be large) absorbed instantly and homogeneously, for example via a chemical reaction with other species present in 
the flow. Therefore, all the predictions of the letter may be checked experimentally as well as numerically. 

Consider advection of passive scalar 6(t;r) by a smooth incompressible velocity field, u(f;r) = <3"(t)r (a(t) is d x d 

. 1 traceless matrix) in a presence of a linear damping and diffusion, 

d t 6 + a*" {ty^O = K& r 6 - a6 + <j>. (1) 

The scalar is forced by random field </>(<; r), which for a sake of simplicity is considered to be Gaussian thus fixed 
completely by ((p(ti; ri)(p (£2; Y2)) = x(l r i — r 2 1 (i 1 —tz), where the function \{r) decays fast enough if r exceeds 
the integral scale L. Although we are sure that many principal results of the paper can be generalized for the case 
of a finite temporal correlations of the strain matrix tr, only the short correlated one allowing a simple derivation is 
considered here. The Gaussian statistics of a is fixed by the pair correlation function ^cr'" 1 (i)a' 3j/ (t')) equal to 

D [(d + 1) 8^8'^ - 8 m 8 vf) - 8^8^} 8{t-t'). (2) 

We start studying the pair correlation function of the scalar field: F{r\2) — (0(t; r±)8(t; r 2 )). Averaging two replicas 
of (0) multiplied by the respective scalar counterparts one gets 



l - d d r r d ( D(d - l)r + — )<>,. 2n 



F(r) = X (r). (3) 
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Consider the case of a step-like pumping function, when x( r ) = P — const at r < L, and zero otherwise. Introduce a 
forced solution of this equation, Fj(r) = Pd{L — r)/[2a\. Two zero modes of the operator from the lhs of (||) should 
be added to Ff(r) to respect continuity of F(r) and its derivative at r = L. Zero mode added at the upper ( r > L) 
interval should vanish at r —> L, while its counterpart added at r < L should be finite at the origin ( r — > 0). If the 
dissipative scale, r d = */ k/ max{Z), a} is small enough, one gets 

r d ^r<L, 

(4) 

r > L > r d ; 




where £± = ±^d 2 /A + 2ad/ [(d — 1) D] — d/2. Finite diffusivity generalization of (Q) (of its zero-mode part) can be 
presented in terms of the hyper-geometric function. Dominant contribution into S^M at r C L stems from a zero 
mode of the operator on the rhs of (||) scaling as . 

Come to the study of the scalar difference stationary PDF, V = (5 (x — S6 r )). Generally, at zero diffusion (k — > 0) 
the stationary limit is perfectly achieved via the direct balance between the pumping <j> and the a-damping. Therefore, 
there is no dissipative anomaly in the case and we may simply write the Fokker-Planck equation 

D(d - iy- d d r r d+1 d r V + ad x (xV) = 0. (5) 

(||) is valid at i < 6l, where 9^ is the amplitude of the scalar field at the integral scale, which is estimated by 
P/max{a, D}. Even without calculation of the PDF itself we may simply get the anomalous exponents for the 
structure functions of all the orders, by means of integration of (^J) against the respective moments of x. The only 
thing left is to pick up a vanishing at r — > solution of the linear ordinary homogeneous differential equation (zero 
mode of the eddy-diffusivity operator) . It gives for the even moments (odd moments are constrained to be zero due 
to pumping isotropy and Gaussianity) 

S2n(r) fr\t*« t Id 2 2adn d 




where the dependence on L and 6l is restored just by dimensionality. (^|) is compatible with (^) and it holds for any 
d , n, and a. Notice that at the fixed value of the damping-to-advection ratio a/D and d going to infinity £2™ goes 
to zero. 

Although the calculation of anomalous exponents was our main goal the PDF itself may be also simply extracted 
out of (||) if additionally the form of the PDF at the integral scale (which in its own term relates to concrete form of 
the pumping) is fixed. For example, the Gaussian PDF at the integral scale corresponds to 

V&r)= l ^(±) d/2 Q(pML/r)), (7) 

dt 
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where a = D(d — l)/a. At L r the PDF shows a change in behavior about x c = 9l [r / L] 1 ^ da \ Q (y; z) is finite at 
the origin, Q (0; z) ~ (r/L) d ^ 2 / In [L/r], and further the alrgebraic in x decay at x <C x c , Q (y; z) — Q (0; z) ~ y d ' a / 4 : 
turns into Q (y; z) ~ y d a / 4 exp [— z 2 / (aln [1/j/])] at x x c . Therefore, (^|) is applicable for all the positive (not 
necessarily integer) moments of \59 r \. All the negative moments are divergent. 

The possibility of two-point consideration explained above is based on the absence of the dissipative anomaly. To 
prove this and also to shed some light on the dynamical origin of anomalous behavior we consider Lagrangian multi- 
point representation of the problem and show how does it lead to the same answer ([}]) for the asymptotic behavior of 
the structure functions. ([!]) is equivalent to 

00 

0(t; r) = / dt' exp [-at'} (f; p(t - t')) 

(°) 
f t p(t) = a(t)p(t)+n(t), p(0)=r, 

where 77(f) is the Langevin noise fixed by, (r] a {t)n^ (t 1 )^ = 2n5(t — t'). Averaging the simultaneous product of 2n 
different replicas of (pi) one gets F\...2 n = {61 ■ ■ ■ #2n)- Once the average of the multi point product is known it is easy 
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to construct the desirable structure function S2n(r) fusing the 2n points. Moreover, one can get S2n from the general 
object with all the points being placed along a straight line. Also we expect, and it will be confirmed below, that 
at a > all the fused averages are finite also in the limit of zero diffusivity. The last observation allows to consider 
the infinite Peclet number Pe = L/ra limit just replacing k by zero. Therefore, for the collinear geometry = n 
and at Pe — * oo direct averaging of (M) with respect to statistics of (f> and a gives 



l---2n 



{l,---,2n} 

E ■ 

{*!,.. .,*2t»} 



n 

\k=l 



dt k e- at « X 



(9) 



where r)(t) = \W(t)n\ and W(t) satisfies dW(t}Jdt = a(t)W(t). The a = version of (§) was calculated in gj] for the 
d = 2 case and generalized for any d > 2 in fl5| via a change of variables and further straightforward transformation 
of the path integral standing for the average over a(t). Average with respect to all the "angular" degree of freedoms 
gives the following effective measure of averaging with respect to the rate of stretching along the n direction |14|,|l5fl 



T>r](t) exp 



-d J dl 
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fi 2 + D 2 (d-1) 2 - 2(ri-l»£> 
4(d- l)D 



Therefore, averaging with respect to the fluctuations of the rj field produces 



d D(d-l)d 
77 ^1 Z tl 



OO tl tn — 1 

Fi-2n = n\ J dti Jdt 2 --- J dt n J drji ■ ■ ■ dr) n exp 

°° 

{l,-,2n} n 

x E II[e 3at4 x(e' K r fca4 , fa4+1 )G(t j _i-t i ;T ?i _i-T h )] 

{ki,-,k 2n }i=l 



(10) 



where £j + i = 77^1 = and G (i; 77) = 



ird 



i{d-X)tD 



exp 



is the Green function of the diffusive kernel. The 



A(d-l)tD 

integrand of (|10|) decays exponentially in time, so the major contribution into the object forms at tj ~ it does not 
depend on any r^ . The first r-dependent contribution stems from n — 1 temporal integrals formed at r ~ 1/a, and 
one at ~ r r ~ In [i/r] /max{a, ZP} (this special integration brings a spatial dependence into the object, therefore 
on a single distance). Generally, there exists a variety of terms with all the possible combinations, like a term with k 
integration formed at r, while n — k ones at r r , and therefore dependent explicitly on 2(n — k) points. However, we 
are looking exclusively for a term dependent on all the 2n points since only such a term of (|10|) contributes <S , 2n(^)- 
And it is really simple to calculate the scaling of the term making use of the temporal scale separation, r r ^> r. 
Indeed, the large time contribution may be extracted out of ( |Io| ) in a saddle-point calculation. Variation of all the 
exponential terms in ( |Io| ) with respect to ti gives a chain of saddle equations. The x functions in the integrand 
of (^) limits the 77 integrations from above by In [L/r]. Therefore, the desirable 2n-points contribution forms at 
ti = y/ dj [4(d — 1)D (2cm + D(d — l)d/4)] In [L/r], and r\i — In [L/r], where in the leading logarithmic order one 
should not distinguish between contributions of different separations . Finally, substituting the saddle-point values 
of ti and rji into ([[o|) we arrive at (||). 

The basic physics of nonzero <^2n (means deviating from the naive balance of pumping and advection) and generally 
anomalous (A2„ 7^ 0) scaling at a > can be stated quite directly. According to (|]) the advection changes scales 
but not amplitude, while the amplitude of injected scalar field decays exponentially from the time of injection at the 
constant rate a. The temporal integrals in (^|) forms at the mean time to reach a scale which is proportional to the 
negative log of the scale. However, the effective spread in the factor by which amplitude has decayed, upon reaching 
a given scale, increases as scale decreases. It is why ^2n > 0. Also there is more room for fluctuations about the mean 
time due to the interference between the exponential decay of the scalar amplitude and fluctuations of the stretching 
rate 77. Thus intermittency increases with decrease of scale size. 

I conclude by couple of general remarks. First of all, the model gives an example of the situation when the dissipative 
anomaly is absent, while the intermittency (anomalous scaling, A2„ ^ 0) takes place. Second, continuous dependance 
of the exponents on the damping rate originates from coincidence of the scaling dimensions (zero in the Batchelor 
case) of the bare eddy diffusivity operator and the damping-dependent correction to it. 

I thank I. Kolokolov, R. Pierrehumbert, B. Shraiman for inspiring discussions. Valuable comments of G. Falkovich, 
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